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ABSTRACT 

We show that a one-component variable-emissivity-spectral-index model (the free-/3 
model) provides more physically motivated estimates of dust temperature at the Galactic po- 
lar caps than one- or two-component fixed-emissivity-spectral-index models (fixed-/3 mod- 
els) for interstellar dust thermal emission at far-infrared and millimeter wavelengths. For the 
comparison we have fit all-sky one-component dust models with fixed or variable emissivity 
spectral index to a new and improved version of the 210-channel dust spectra from the COBE- 
FIRAS, the 100 - 240 /im maps from the COBE-DIRBE, and the 94 GHz dust map from the 
WMAP. The best model, the free-/3 model, is well constrained by data at 60 — 3000 GHz over 
86 per cent of the total sky area. It predicts dust temperature (Tdust) to be 13.7—22.7 (±1.3) K, 
the emissivity spectral index (/3) to be 1.2 — 3.1 (±0.3), and the optical depth (r) to range 
0.6 - 46 X 10~^ at = 900 GHz (Aq = 333/im) with a 23 per cent uncertainty. Using 
these estimates, we present all-sky evidence for an inverse correlation between the emissiv- 
ity spectral index and dust temperature, which fits the relation (3 = 1/{S -\- uo ■ Tdust) with 
S = -0.510 ± 0.011 and uj = 0.059 ± 0.001. This best model will be useful to cosmic 
microwave background experiments for removing foreground dust contamination and it can 
serve as an all- sky extended-frequency reference for future higher resolution dust models. 

Key words: dust, extinction - infrared: ISM - submillimetre: ISM - Galaxy: general - meth- 
ods: data analysis - technique: spectroscopic. 
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1 INTRODUCTION 

An accurate model of thermal dust emission at the far-infrared and 
millimetre wavelengths is important for cosmic microwave back- 
ground (CMB) anisotropy studies because it helps to remove one of 
the three major diffuse foreground contaminants. In the last decade, 
experime nts such as the Wilk inson Microwave Anisotropy Probe 
(WMAP, iBennett et al.ll2003ah have precisely measured the angu- 
lar variations in CMB signal in order to understand the global ge- 
ometry and expansion of the universe. However, studying varia- 
tions that are 10~^ the strength of the principal signal is difficult, 
and the removal of contaminating signals in the data needs to be 
done accurately. For t hese experiments, a dust templa t e, such as one 
extrapolated from the iFinkbeiner. Davis & Schlegell (Il999h (FDS) 
study of interstellar dust in the far-infrared has been u sed to remove 
thermal dust contribution from sky measure ments jBennett et al.l 
l2003bl : lHinshaw et al.ll2007l : lGold et al.ll2009h . 

Among major efforts to derive an all- sky dust model from 
observational data, iReach et aP (Il995h use dust spectra derived 
from measurements of the Far Infrared Absolute Spectrophotome- 



E-mail: zliangl@jhu.edu 



ter (FIRAS, IfIRAS Explanatorv Supplement! (Il997l)) instrumen t 
on board the Cosmic Background Explorer ( COBE. iMatheJl 1 982h 
satellite to constrain dust models with emissivity proportional to 
ly^. They found that dust emission is best described by a three- 
component dust model: a warm (16 — 21 K) and a cold (4 — 7 K) 
component that are present everywhere in the sky, and an interme- 
diate temperature (10 — 14 K) component that exists only at the 
Inner Galaxy. In 1996, Boulanger et al. independently derived an- 
other set of dust spectra using the FIRAS measurements and fit to 
it a one-component emissivity dust model. They found that the 
average spectrum of dust associated with HI gas had an averaged 
temperature of 17.5 ± 0.2 K. In 1998, Lagache et al. used DIRBE 
bands at 100, 140 and 240 /im to decompose FIRAS spectra at 
\b\ > 10° into a cirrus and a cold component. For 61 per cent of 
the sky where the cold emission was negligible, they found that the 
cirrus had a mean temperature of 17.5 K with a dispersion of 2.5 K. 
For the 3.4 per cent sky where both cirrus and cold components 
were present, the two components were both assumed to follow a 
emissivity law, with the cirrus component found to have a tem- 
perature of 17.8 ± 1.2 K, and the cold component with a temper- 
ature of 15 ± .8 K. A widely used du st model in CMB studies 
was obtained bv IFinkbeiner et al.l (Il999h . Their best model (Model 



©2011 RAS 



2 Z Liang, D. J. Fixsen and B. Gold 



#8) to the FIR AS dust spectra consisted of two components: a cold 
component following a zy^^^ emissivity law with temperature at 
7.7 — 13.1 K, and a warm component following a zy^^° emissivity 
law with temp erature at 13.6 — 21.2 K. A d ecade later, the Planck 
Collaboration dPlanck Collaboration! [20 1 1 ah used the Planck-UFl 
(350 /im-2mm) and IRAS 100 /am data to derive all- sky dust tem- 
perature and optical depth maps using a one-component model with 
emissivity proportional to z/^'^. They found that the median tem- 
perature of the sky at 10° abo ve and below the Galactic plane was 
17.7 K, see also iLiand (1201 ih . 

The above results highlight the diverse findings in the study of 
thermal dust emission at far-infrared and millimetre wavelengths. 
They show that the derived dust properties depend as much on the 
fitting method and the functional form of the model as on the data. 
With added new and more sensitive data from the WMAP, we now 
can constrain model parameters with much greater accuracy. 

A second reason for our work is to present fit results from 
applying variable spatial averaging to increase signal-to-noise of 
the spectra. Since low intensity measurements often come with 
large uncertainties, when such data are used directly to constraint a 
model, results are highly uncertain parameters. At times this prob- 
lem is treated with averaging data within a predefined sky region. 
This approach has the disadvantage of using a presupposed dust 
distribution in the derivation of a solution while figuring out the 
distribution is part of the research question. Here we make no as- 
sumption of the dust distribution but instead use the signal-to-noise 
of the data to determine the amount of spatial averaging needed for 
the data. The results are higher spatial resolution for regions with 
good signal-to-noise and less averaging for the original data set. 

A third motivatio n for our work is to understand whether 
dust optical properties jPraine & Led 1 19841) are the same at far- 
infrared and millimetre wavelengths from the perspective of em- 
pirical model fitting. That far-infrared dust emissivity follows a 
power law has been widely accepted (see list above), yet the va- 
lidity of such an extrapolation has not been proven by theory, lab- 
oratory experiment, or empi rical model fit ting. In fact, reports of 
labora t ory measureni e nts by lAdadze et al. (ll996l):lMennella et al.l 



labora t ory measureni e nts by lAglaaze et al.l (li^S^ol):IMennella et al.l 
(Il998h : iBoudet et all (120051) : ICoupeaud et all ( 1201 ll) have shown 



that emissivity of amorphous silicate and carbon grains differed 
from a power law and had a significant temperature dependence. 
This inconsistency between laboratory measurements and model- 
ing of astronomical observations means that our understanding of 
dust emission in the far-infrared and millimetre is incomplete. In 
this work, we attempt to fill this gap by first deriving best-fitting 
dust models with emissivity spectral index fixed at different val- 
ues and as a variable, and then comparing the quality-of-fit of these 
models. Based on our findings, we argue that dust emissivity differs 
in the far-infrared and millimetre from the optical. 

Finally, regarding the many empirical models we now know, 
e.g. the ones listed above, one cannot help but ask: How do we 
test the validity of these models? Beside having good constraints 
on model parameters, are there physically motivated tests we can 
use to verify predictions of these models? Here, we propose one: 
to compare dust temperature distribution with the distribution of 
dust heating source at the Galactic poles. We conduct an indepen- 
dent and comprehensive test on all-sky one-component models, and 
show that all but the one-component free-/3 model fail this test. 

The structure of our manuscript is as follows. In Section[2l we 
review observations by the COBE satellite's DIRBE and FIRAS 
experiments and the WMAP satellite that are used in our model 
construction. In Section [3] we detail procedures taken to deduce a 
new set of FIRAS dust spectra and to unify calibrations of the data 



sets. In Section|4l we present results and analysis from fitting one- 
component dust models with fixed and variable emissivity spectral 
index to the data. In Section^ we compare temperature predictions 
of the free-/3 model with those of the fixed- /3 models and show 
that only the free-/3 model gives physically motivated predictions 
of dust temperature at the Galactic polar caps. We also discuss the 
implications of the free-/3 model on the inverse correlation between 
emissivity spectral index and dust temperature. Conclusions along 
with suggestions for how to use our results are presented in Section 



2 OBSERVATIONS 
2.1 COBE DIRBE 

The DIRBE instrument was a cryogenically cooled 10-band ab- 
solute photometer designed to measure the spectral and angular 
distribution of the diffuse infrared background. It had a 0?7 beam 
and covered the wavelength range from 1.25 to 240 /xm. During its 
lifetime, the DIRBE achieved a sensitivity of 10~^ W m~^ sr~^ 
at most wavelengths (iBoggess et aDll992l : Isilverberg et al.llT993l : 
IdIRBE Explanatorv Supplemenj|l99"8h . 

We use the 1997 "Pass 3b" Zodi-Subtracted Mission Aver- 
age (ZSMA) Maps at bands 100, 140, and 240 /xm. These maps 
measure the Galactic and extragalactic diffuse infrared emission 
and have been calibrated to remove zodiacal light (zodi). They are 
available at the Legacy Archive for Microwave Background Data 
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Analysis (LAMBDA 



2.2 COBE FIRAS 

The FIRAS instrument was a polarizing Michelson interfer- 
ometer designed to precisely measure the difference between 
the CMB and a blackbody spectrum. The FIRAS had a 7° 
beam and covered the frequency range from 1 — 97 cm~^ at 
0.45 cm~ resolution (Boggess et alj 1 19921 : iFixsen et al.l Il994al : 
IfIRAS Expla natorv Supplement 1993. 

We derive a new set of dust spectral maps from the Destriped 
Sky Spectra of the Pass 4 final data release using procedures de- 
scribed in Section [3] The 210 6063-pixel maps comprise the main 
body of spectral information in our model fitting. 

Six types of uncertainties have been characterized by the FI- 
RAS Team The instruct i ons on h ow to treat each of them are doc- 
umented in 'Fixsen et al.l Jl994bl) : rFIRAS Explanatorv Supplement 
(1997); Mather et al. (1993). Since we want to build models that 
respond to both spectral and spatial variations of dust emission, 
our analysis includes all six FIRAS uncertainties: detector noise 
(D), emissivity gain uncertainties (PEP), bolometer parameter gain 
uncertainties (JCJ), internal calibrator temperature errors (PUP), 
absolute te mperature errors (PTP), and destriper e rrors (p). Sec- 
tion 7.10 of IfIRAS Explanatorv Supplement! (Il997l) provides very 
helpful instructions on how to assemble the covariance matrix. For 
example, the D and /3 matrices vary only among pixels, while the 
PEP, JCJ, PUP and PTP matrices differ for different frequencie s. In- 
tereste d readers are referred to the iFIRAS Explanatorv Supplement! 
(!l997!) for details. 



^ The LAMBDA Web site is http://www.lambda.gsfc.nasa.gov/ 
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Table 1. Spectral Coverage of DIRBE, FIRAS & WMAP 
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2.3 WMAP 

The Wilkinson Microwave Anisotropy Probe was designed to de- 
termine the geometry, content and evolution of the universe by 
measuring temperature anisotropy of the CMB radiation. It con- 
sisted of two back-to-back offset Gregorian telescopes and used 
20 high electron mobility transistor (HEMT) based differential 
radiometers to measure the brightness difference between two 
lines of sight that were 141° apart. At five frequency bands: 23, 
33, 41, 61, and 94 GHz, the WMAP made full sky measure- 
ments, which were analyzed by the data processing pipeline and 
formed 13^ FWHM HEALPi^Q pixelization maps. The spin mo- 
tion of the observatory and its scanning strategy symmetrized 
the WMAP beams. Beam sizes were estimated using square- 
root of the beam solid angle. In the order of increa sing frequen- 
cies, t hey were: 0.88, 0.66, 0.51, 0.35, and 0.22° (j arosik et alj 

Jarosik et al. 2003; 



20031: IPage et"aD l2003l: iHinshaw et all [20031; 
WMAPFive-Year Explanatorv Supplement! 



2009l : lHilletal.ll2009l) . 



20081: IHinshaw et all 



We use the dust temperature map (at 94 GHz) derived from 
the "base modeFin WMA P's Five- Year foreground modeling anal- 
ysis by Gol d et all (l2009h . In the same study. Gold et al. used 
different models to account for the diffuse foreground emission 
at different WMAP bands, with nonthermal synchrotron, thermal 
bremsstrahlung, and thermal dust as the standard components and 
tested the possible existence of steepening synchrotron and/or spin- 
ning dust. Their likelihood analysis has shown that the basic model 
with just three main foreground components would be sufficient to 
remove foregrounds from sky maps at high Galactic latitudes. 



3 DATA PREPARATION 

Since our objective has been to model both spectral and spatial vari- 
ations of dust emission, we unified different hardware constraints 
and calibration standards to ensure that different data sets are com- 
pared on an equal footing. In Section lTTI we explain the deduction 
of FIRAS dust spectra. In Sections [T2l - 1 3. 7 1 we explain treatments 
on the DIRBE and WMAP data to make sure that they have the 
same physical attributes as the FIRAS dust maps. Among proce- 
dures presented in this section, those pertaining to beam differ- 
ences, map projections, spatial resolutions, DIRBE-FIRAS abso- 
lute calibrations, and temperature-flux conversion are applied to 
both signal and noise maps. Procedures on zodi zero-point correc- 
tions and FIRAS systematic errors are applied to noise maps only. 
The uncertainties used in our analysis are the quadratic sums of the 
uncertainties of data sets and models described in this section. More 
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Figure 1. Examples of the new FIRAS dust spectra. Plotted in red are the 
new dust spectra and their error estimates; in purple are the dust spectra 
derived by the FIRAS team; in dark green is the total sky intensity measured 
by FIRAS; in blue is the CMB monopole; and in orange is the CMB dipole. 



specifically, when we are given uncertainties of the parameters, we 
use 

2 



3 ^ ^ 3 k^3 



to generate uncertainties of intensity predictions. 



3.1 Deducing FIRAS dust spectral maps 

The FIRAS dust spectra on LAMBDA exhibit a "jump" between 
the low- and high-band data due to an inconsistent CMB monopole 
temperature subtraction. In addition, the Pass 4 data were released 
with an early calibration, making it necessary that we derive a new 
set of dust spectra. In the following we describe procedures to sub- 
tract from the Destriped Sky Spectra a blackbody spectrum for the 
CMB, a dipole of the Earth's motion with respect to the CMB, a 
zodi model, and contribution from the cosmic infrared background 
(CIB). Examples of the new dust spectra are plotted in red in Fig. 
[11 For comparison, corresponding dust spectra provided by the FI- 
RAS Team are plotted in purple. 



3.1.1 CMB monopole and dipole 



The CMB temperatur e has be e n extensively t r eated i n lMather et al 



19901) : iFixsen et all ( Il994bh : |^^^^^J^J^_E^^^ 
199^ : iMather et al.l (Il999l) : iFixsen & Dw"e3 (l2002h : Fixsed 
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(l2009h . The appropriate correction for the Pass 4 data set is a 2.7278 

black body spectrum. 

A WMAP-determined dipole (iHinshaw et al.ll2009l) has been 
removed from the Destriped Spectra. Specifically, Tdipoie = 
3.355 mK and (/, 6) = (263? 99, 48?26). Higher order variations in 
the CMB temperature were ignored because they were insignificant 
for this study. 

3.1.2 Zodi 

Zodi is the thermal emission and scattered light from interplan- 
etary dust in our solar system. iKelsall et"al I (Il998h derived a 
time-dependent parametric model for its emission, and the FI- 
RAS Tea m extended those results to the entire FIRAS frequency 
cover age (IfIRAS Explanatorv Supplement! 1 19971 : iFixsen & Dwekl 
l2002h . Derivation of the zodi model for FIRAS hinged on the fact 
that FIRAS measurements overlap with DIRBE bands at 140 and 
240 /im. Therefore, the DIRBE model predictions for these two 
bands were coadded and fitted with a power law emissivity model 
and extrapolated to the frequency coverage of FIRAS. The zodi 
model used here is among FIRAS data pro ducts on LAMBDA. Fo r 
further details of the model derivation, see lFixsen & Dwe3 (l2002h . 



3.1.3 Emission lines 

The FIRAS detected 1 8 molecular and atomic lines emitted by in- 
terstellar gas. Since the FIRAS frequency resolution is much larger 
than the width of each of these lines, each line profile is effectively 
FIRAS 's instrument response to a delta function. Among these 18 
detected emission lines, not all of them have a discernible presence 
over the full sky. Most notable are the [C ii] and [N ii] lines, which 
exhibit a distinct gradient of intensity from the centre of the Galaxy 
to higher latitudes. Other emission lines, though detected, are weak 
in most of the sky except at the Inner Galaxy. By Inner Galaxy we 
mean the inner Galactic disk about half the distance to the edge of 
the Galaxy. 

Since the derivation of FIRAS line intensity maps on 
LAMBDA used their Galactic dust spectra, the FIRAS line inten- 
sity maps cannot be used here to remove emission line contribution. 
Instead, intensities of [C ii] and [N ii] emission were fit as parts of 
the overall model. 



3.1.4 Cosmic infrared background 

We removed the isotropic CIB signal from sky spectra by using 
results from three studies: For DIR BE measurements a t 140 and 
240 /im, we used estimates given bv lHauser et al at 15 and 

13 nW m"^ sr~ \ For the DIRBE band at 100 /im, we use the 
estimate given in Finkb einer et "Zl boooh at 25 nW m ^ sr ^ , see 
also Lagache et al. (200^. To re move the CIB froni FIRAS sky 
spectra, we used the CIB model in lFixsen et al.l (Il998h . 



3.2 Beam difference 

The three instruments that produced the data used in this study 
had different beam patterns. For example, the FIRAS used a quasi- 
optical m ultimode horn antenna to collect radiation from a 7° field 
of view (M ather et al.|[l986i) . The horn was designed in a trum- 
pet bell shape to reduce response to off-axis radiation. As a result, 
when the beam profile was measured on the ground and in flight, 
it was found to have very low sidelobes over the two decades of 



frequency measured by FIRAS. The central portion of the beam 
{0 < 3? 5) was approximated by a top hat since any slight azimuthal 
asymmetry should have been smoothed out by th e rotation of the 
instrument along its own axis jMather et al.lll993h during its oper- 
ation. 

On the other hand, DIRBE was built with a goal to reject 
stray light to measure the absolute spectrum and angular distribu- 
tion of the CIB. This goal was met by using a series of optical 
elements and baffle protections, among which the last field stop 
set the 0? 7 x 0? 7 instantaneous field of view for all spectral b ands 
dSilverberg et al.ll 19931 : IdIRBE Explanatorv SuDplementlll998h . To 
construct the ZSMA maps, the DIRBE Team calculated the zo dia- 
cal light intensity using the IPD model by Ke lsall etaP ((1998) and 
subtracted it off from each weekly measurement. The remaining 
signal was averaged over time. In this way, the ZSMA maps pre- 
served the original 0? 7 x 0? 7 angular resolution of the sky obser- 
vation. 

The dust map from the WMAP was one of the products de- 
rived from Marko v chain Monte C arlo fitting of temperature and 
polarization data (iGold et al.11260^ . Since their analysis used the 
band- averaged maps that were smoothed by a 1 ° Gaussian beam, 
the dust map had the same angular resolution. 

Since the FIRAS beam had the lowest common angular reso- 
lution achievable among all three data sets, we convolved the higher 
resolution DIRBE and WMAP maps w ith the FIRAS beam t o make 
them all 7° maps. Additionally, since iFixsen et al.l (Il997bh found 
that the FIRAS beam was elongated in the scan direction by 2? 4, 
we matched that pattern in the degraded DIRBE and WMAP maps 
by convolving those data with an effective FIRAS beam. 



3.3 Map projection and spatial resolution 

Both DIRBE and FIRAS maps are organized in CO BE quadri- 
lateralized spherical cube format (quad-cube, tChan^ O'Neilll 
1971 lO'Neill & Laubschejll976L IWhite & Stemwedellll99l and 



Calabretta & Greisenll2002h . While the DIRBE maps are in quad- 
cube resolution level 9 (res9, 19^43 per pixel), the FIRAS maps 
are in quad-cube resolution level 6 (res6, 2? 59 per pixel). Different 
from the F IRAS or the DIRBE maps, the WMAP dust map is in 
HEALPix (iGorski. Hivon & Wandeltl[T999llG6rski et al.l l2005. and 



ICalabretta & Roukemal l2007h resolution level 6 (res6, 54™ 97 per 
pixel). One way to reconcile these different formats and spatial res- 
olutions is to carry out analysis in COBE quad-cube res6. With this 
decision we hope to retain maximum amount of information con- 
tained in the original data sets and to achieve the highest common 
resolution possible. 

We re-binned DIRBE maps to res6. The WMAP dust map was 
first converted into a quad-cube res9 map and then re-binned to 
res6. During WMAP's dust map conversion, we checked to make 
sure that no excessive artificial noise was introduced to the final 
map: By comparing the original HEALPix-projection map with the 
re-binned quad- cube-projection map, we found that 98.6 per cent 
of the 49,152-coordinate pairs sampled gave no difference be- 
tween the quad-cube and the HEALPix values. When there was 
a difference, the maximum was 0.0059 mK, which amounted to a 
0. 1 1 per cent noise increase for the original HEALPix map. 



3.4 Gradient correction 

Since the production of FIRAS dust maps involved coadding in- 
terferograms, the value reported for each pixel is generally defined 
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for a location within the pixel other than its defined centre. This 
positional difference requires additional correction to prepare the 
quad-cube res6 DIRBE and WMAP maps. Details of this technique 
are described in lFixsen et al.l (Il997bh . In summary, we fit a second- 
degree surface function to the intensity and location information of 
a pixel and its immediate neighbors in one of the converted maps. 
We then use this function to predict emission at the FIRAS mean 
position for that particular pixel. Overall, we measure a 5 per cent 
rms correction to each of the DIRBE maps at 100, 140 and 240 /xm 
and to the WMAP dust map. 



3.5 Color correction 

In accordance with the IRAS convention 

dlRAS Explanatorv Supplement! Il988h . DIRBE photometric 
measurements were reported in MJy sr~^ at nominal wave- 
lengths, assuming the source spectrum to he v • 1^ = constant. 
Since each DIRBE band has a much wider bandwidth than a 
FIRAS channel, spectral shape could have changed enough that 
at the nominal wavelength the real intensity was significantly 
different from the normalized intensity. As a result, we include 
color corrections in the overall model, i.e., model predictions 
are compared with DIRBE measurements using the relation 
lu, model = lu, DIRBE /i^, whcrc K IS the color correction factor 
defined as 



K : 



J (z/o/i^)quoted ' Ru 



(1) 



In this equation, J (/iy//i/o) actual is the specific intensity of the 
sky normalized to the intensity at frequency and is DIRBE 
relative sys tem response at frequency u. The values of Ri, are doc- 
umented in lDIRBE Explanatorv Supplement! (Il998h Section 5.5. 



3.6 DIRBE uncertainties 

The DIRBE photometric system was maintained to ^ 1 per cent 
accuracy by monitoring the internal stimulator during 10 months 
of cryogenic operation and observing the bright stable celestial 
sources during normal sky scans. It was absolutely calibrated 
against Sirius, NGC7027, Uranus, and Jupiter. 

Among different types of uncertainties identified by the 
DIRBE Team, those relevant to this work are standard devia- 
tions of intensity maps, de tector gain and offset uncertainties, and 
zodi model uncertainties jHauser et al ] |l998l : iKelsall et al.lll998l : 
lArendt et aDE998). For bands 8-10, respectively, we use detector 
gains estimates at 0.135, 0.106 and 0.116 nW m~^ sr~^, detector 
offsets at 0.81, 5 and 2 nW m~^ sr~^, and zodi mo del uncertainties 
at 6, 2.3 and 0.5 nW m"^ sr"^ (lArendt et al.lll998h . The process of 
re-binning the high resolution DIRBE maps into FIRAS resolution 
affects the standard deviations of the original maps only. The final 
uncertainty is the quadrature sum of the individual noise compo- 
nents. 



3.7 Temperature-intensity conversion 

Foreground maps of the WMAP production are reported in antenna 
temperature, Ta, in mK. On the other hand, maps produced by the 
DIRBE and the FIRAS Teams are reported in spectral intensity, /jy, 
in MJy sr~^. In the following analysis, the WMAP dust map is 
converted into flux density values following 1^ — 2 (y/cf kTA, 
where v is the effective frequency (93.5 GHz) of the dust map 



dOold et al.|[2009l : Ijarosik et al.ll2003h and k is Boltzmann's con- 
stant. 



4 RESULTS AND ANALYSIS 

4.1 Overview of model fitting strategy 

The thermal emission of a dust grain can be modeled by a modified 
blackbody function: 

^dust(^') = r Ely ' Bjy{Tdust) , 

where Bu{Tdust) is the blackbody spectrum at temperature Tdust, 
= (y /vqY is the emissivity with spectral index /3, and r is the 
optical depth normalized to frequency z^o = 900 GHz. 

In addition to measuring thermal dust emission, the prepared 
FIRAS spectra retain contributions from [C ll] and [N ii] emission 
due to the lack of precise all- sky templates. As a result, the two 
emission lines are modeled as: 



and 



(l/) = [Cii] intensity /[C„](^) ' 



^[Nii](2^) = [N II] intensity /[N„](^)5 

where f(v) is the synthetic line profile determined by FIRAS re- 
sponse to a delta-function signal. Together, the full model has the 
form 



-''total — I dust + -^[Cn] + -^[Nii] • 



(2) 



Each full model is fit to the data by minimizing a three-part 



X with each part corresponding to one of the three data sets: 

2 _ 2 I 2 ,2 

X — XdIRBE + XfIRAS + XWMAP 5 



(3) 



where 



Xinstrument 



y^(Jobs — Imdlji (M (/obs — Imd\)j , (4) 



Here, Jobs is the observed spectral intensity, Imdi is the model pre- 
diction, and M is the covariance matrix of the respective data set. 

In the following sections, we fit one-component dust models 
to spectra of fixed (Section |4l2l ) and different (Section |4]3]) size sky 
regions. In the former case, spectra retain the 7° angular size of FI- 
RAS pixels; in the latter case, the 7° spectra are averaged by various 
amounts to increase signal-to-noise of the final spectra. Chi- square 
per degree of freedom values, Xdof » ^re used to assess the quality- 
of-fit of a model to each spectrum. In particular, the one-component 
fixed- /3 models have 214 — 4 = 210 degrees of freedom while the 
free-/3 model has 209 degrees of freedom. Adopting a 10 per cent 
probability cut-off, it corresponds to a Xdof ^ 1.13 for both fixed- 
and free-/3 models. 



4.2 Fitting spectra of 7° sky regions 

We fit to the spectrum at each 7° pixel one-component models with 
fixed 13 in the range 1.4 — 2.6 at 0.1 increment. The fits have ac- 
ceptable values (Xdof ^ 1-13) over most of the sky except at the 
Galactic plane. 



4.2.1 Quality -of -fit of models 

As an example. Fig. |2] presents Xdof ^s a function of Galactic lat- 
itude for the P = 2.0 model. It shows that most fits at |6| > 10° 
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Figure 2. x^of Galactic latitude. The value of x^of obtained from 
fitting one-component (3 = 2.0 model to 7° spectra covering 98.68 per cent 
area of the full sky where FIR AS data are available. 



250 




Figure 3. x^of distributions of the 7° fits using one-component (3 = 2.0 
model. The red distribution includes only pixels at Galactic latitudes |6| > 
10°; the blue distribution includes all 6063 pixels at all Galactic latitudes. 
Best-fitting parameters of the Gaussians are printed in respective colors. 



have Xdof ~ 1 that fits at |6| < 10° have a Xdof > 1- Fig- 
[3] compares the distribution of Xdof \b\ > 10° with the distri- 
bution of Xdof for the entire sky. That both distributions are well 
approximated by a Gaussian indicates that the fits don't have a sig- 
nificant systematic bias. The widths of the distributions are as ex- 
pected ( 0.1) for a distribution of random data with 210 degrees of 
freedom. That the distributions center at 0.93 means that statistical 
errors of the data are slightly overestimated by ^ 7%, and that the 
Xdof cutoff is really at ^ 1.21 with a probability of < 10%. Be- 
cause the uncertainties of the FIRAS data include some systematic 
effects, we do not feel at liberty to reduce the uncertainty. Based 
on the Xdof ^113 cut, the model is a good fit to the data for over 
87 per cent of the full sky area and is rejected by the data at the 
Galactic plane. 

Readers interested in the best-fitting parameters (Tdust, r, 
[Cll] and [Nil] intensities), their uncertainties and c orrelat i ons fo r 
each of the aforementioned P models are referred to lLiangI (1201 ih . 



4.2.2 Dependence of best-fitting parameters on the 
signal-to-noise of data 

Figs. H] and [Upresent Xdof' ^dust and r of the best-fitting {5 models 
for two spectra: one has high signal-to-noise (HSN) and is at a low 
Galactic latitude; the other one has low signal-to-noise (LSN) and 
is at a high Galactic latitude. For both cases, plots of Xdof ^s. (3 
show that models with a wide range of different p values can fit the 
data well. In the HSN case, a curve fit to Xdof as a function of /3 is 
a concave up parabola, with the minimum Xdof — 0-89 at /3 1.8. 
The difference between Xdof at /3 = 1.8 and that at /3 = 2.0 is 
AXdof = 0.01. At 210 de grees of freedom, this means a Ax^ of 
^ 2.1, a 2-sigma difference. In the LSN case, the best-fitting curve 
to Xdof ^s. /3 is a much flatter parabola over 1.4 ^ /3 ^ 2.3 with 
the minimum Xdof = 0-80 at /3 = 1.6. The difference between 
X^^f at /3 = 1.8 and /3 = 2.0 is - 0.002. 

Although models with different (5 have only a small difference 
in Xdof' the best-fitting Tdust and r are different significantly in the 
HSN case: Pd {5 — 2.0, the best-fitting dust temperature is Tdust = 
17.5 ± 0.26 K, compared to Tdust = 18.5 =b 0.28 K at /3 = 1.8. 
This difference in temperature is larger than the sum of their errors. 
Similarly, the difference in r of the two (3 models is larger than the 
sum of their errors. On the contrary, in the LSN case, the difference 
in the best- fitting Tdust and r of p = 1.8 and 13 = 2.0 models are 
well within the uncertainties of the respective parameters. 

These results demonstrate the sensitivity of the fits to measure- 
ment errors. The existence of measurement noise inevitably causes 
a high degree of degeneracy between emissivity spectral index and 
dust temperature in the fits. While it is difficult to break this degen- 
eracy, high signal-to-noise data help. Fitting data with high signal- 
to-noise results in well constrained parameters, which means that 
the choice of an (3 model can cause statistically significant differ- 
ences in the predictions of these parameters. On the other hand, fit- 
ting low signal-to-noise spectra results in small difference in Xdof 
and large error bars of the best-fitting parameters. The results can- 
not be used to differentiate models with different fixed values of 
13. This is demonstrated in Fig. [6l which shows that the 68 and 
95 per cent confidence contours of a HSN fit enclose much smaller 
regions in the T-/3 space than those of a LSN fit. 

Figs. [7] shows sky maps of (3 and Tdust that correspond to 
the minimum-Xdof model among all models tested at each pixel. 
Both maps show greater consistency in value at low latitudes and 
more fluctuations around the Galactic poles. That consistent values 
appear in the region surrounding the Inner Galaxy is reasonable 
because star formation as the most important heat source for dust 
takes place in the Galactic disk and at the bulge. That large fluctua- 
tions appear near the poles, on the other hand, has to do with lower 
signal-to-noise data in these regions compared to those measured 
at lower latitudes. This happens because few dust grains exist at 
high latitudes, and they do not emit as strongly as those close to the 
Galactic disk. 

Since low signal-to-noise data cannot give adequate constraint 
to model parameters and exacerbates the degeneracy between dust 
temperature and spectral index, in order to construct the best dust 
model, we need to find ways to increase the signal-to-noise of the 
data. 



4.3 Fitting averaged spectra of different-size sky regions 

Taking average of the high-latitude spectra based on latitudinal or 
longitudinal divisions of the sky can tighten the constraint on model 
parameters. However, such divisions are based on our expectations 
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Figure 4. Xdof ^dust and r as a function of (5. Each data point on 
the Tdust and r plots is the best-fitting value of the corresponding one- 
component fixed-/3 model to the 7° spectrum measured in the direction 
/ = 63? 78 and h = —11? 53. This set of plots serves as an example of 
high signal-to-noise fits. The green curve fits the best-fitting values as a 
function of /3. 



Figure 5. x^of ' ^dust and r as a function of /5. Each data point on 
the Tdust and r plots is the best-fitting value of the corresponding one- 
component fixed-/? model to the 7° spectrum measured in the direction 
I = 254? 32 and b = 65? 08. This set of plots serves as an example of low 
signal-to-noise fits. The green curve fits the best-fitting values as a function 
of (3. 



of the distribution of Galactic dust. They are not optimal because 
our kn owledge of th e dust distribution is incomplete. In our experi- 
ments (lLiandl20Tll) , of regional fits are much higher than of 
fits to the individual 7° spectra that comprise the regional averages. 
Since larger sky regions include different types of dust emission 
spectra, the steep increase in x^ value means that the averaging has 
achieved a sufficient signal-to-noise ratio and that spectral variation 
has become statistically important. In order to preserve information 
on spectral variation in the model, the amount of spatial averaging 
needs to reflect the signal-to-noise of the data. 

One way to do this is to base the amount of spectral averag- 
ing on signal-to-noise of the averaged spectrum. Starting with the 
base level, where a pixel's own spectrum is used to fit a model, if 
the fit does not give well constrained parameters due to inadequate 
signal-to-noise, the procedure goes on to fit the average of the orig- 
inal spectrum and its eight immediate neighbors. This process of 
involving more of the adjacent spectra to form a new average goes 
on until the derived parameters are sufficiently constrained. In this 
way, results from fits done at the base level have a spatial resolu- 
tion of 6.71 D^. At the next level, results have a spatial resolution 
of 60.37 D^, and so on. 



4.3.1 Tdust/^Tdust constraint on fixed- /3 models 

All-sky one-component fixed-/? models with /3 in the range 1.4 — 
2.6 and lower limit of Tdust / ^^dust at 5, 10, 20 and 40 are obtained 
separately. In general, one-component fixed-/? models with differ- 
ent lower limits on the Tdust /^^^dust values can fit most spectra ex- 
cept those at the Galactic plane. A more restricted lower limit on the 
Tdust /^Tdust values requires a greater amount of spatial averaging. 
As dust spectra with refined differences in shape are binned to- 
gether and the uncertainties of the averaged spectrum are reduced, 
we subject the comparison between the averaged spectrum and the 
modified blackbody function to an increasingly restricted standard. 
The result is a steep increase in x^ value. The lower plot of Fig. [8] 
demonstrates this relation by plotting Xdof as a function of Galac- 
tic latitude for the case of /? = 2.0. Note that Xdof of high-latitude 
(|6| > 60°) fits move from 0.7 - 1.0 to 0.8 - 1.3 as lower limit on 
Tdust /^Tdust starts with none and increases to 40. 

Histograms of the all- sky collections of Xdof for different lim- 
its on Tdust /^Tdust are presented in the upper plot of Fig. [8] For 
the case of Tdust /^Tdust ^ 5, 10 and 20, the three histograms have 
the shape of a Gaussian, peak at 0.93, 0.94 and 0.96 respectively, 
and all have a width of 0.10. At Tdust /^Tdust ^ 40 the histogram 
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Figure 6. 68 and 95 per cent probability contours in the T^^st-f^ space 
for a high signal-to-noise spectrum (upper plot, measured in the direction 
/ = 63? 78 and h = —11? 53) and a low signal-to-noise spectrum (lower 
plot, / = 254? 32 and h = 65? 08). In both plots, the white cross repre- 
sents location of the minimum x^; the blue area is the 68 per cent con- 
fidence region; and the green area encloses the 95 per cent confidence 
region. These plots demonstrate the effect of measurement noise on the 
degeneracy between /3 and T^ust in spectral model fitting. For the high 
signal-to-noise spectrum (upper plot), the 68 per cent confidence region is 
at 1.7 < /5 < 1.9 and 17.9 K < Tdust < 19.5 K; for the low signal-to- 
noise spectrum (lower plot), the 68 per cent confidence region has a much 
wider extent, at 1.2 < (3 < 2.2 and 18.5 K < Tdust < 24.7 K. 




16.0 K 26.0 



Figure 7. Sky maps of /3 and T^ust- We fit all 6000 spectra with one- 
component fixed-/5 models, where 1.4 ^ (3 ^ 2.3 at 0.1 increment, identify 
the minimum-x^ model for each spectrum, and use those models to con- 
struct these maps. In Galactic coordinates Mollweide projection, the centre 
of each map is the Galactic centre. The upper and lower ends of the mi- 
nor axis are +90° and —90° latitudes respectively, and the left and right 
ends of the major axis represent +180° and —180° longitudes respectively. 
Pixels that correspond to fits with a less than 10 per cent probability are 
masked in white. The group of black pixels that slant from the centre of the 
upper left quadrant (North Ecliptic Pole, NEP) to the centre of the lower 
right quadrant (South Ecliptic Pole, SEP) are positions where FIRAS did 
not provide data. 



peaks at LOO, has a width of 0.12, and a thick tail in the range 
1-2 < Xdof < 1-4. This means that at the level of Tdust /^Tdust ^ 
40, it is no longer adequate to use a modified blackbody function to 
describe variations in the dust spectral shape. 

Imposing a more stringent limit on Tdust /^Tiust leads to a 
variety of spatial resolutions in each all- sky collection of fits. A 
more stringent requirement on Tdust / ^Tdust means lower spatial 
resolutions for fits at high latitudes. For the one-component /3 = 
2.0 model with Tdust/^Tdust ^ 5, only 0.59 per cent of the total 
sky area require fits with a 60.37 resolution instead of the default 
6.71 D^. For a 10 per cent constraint on Tdust, 0.47 per cent area of 
the full sky require fits to be at 167.70 resolution, 7.24 per cent 
at 60.37 resolution, and the rest at 6.71 resolution. A balance 
between having adequate constraint on parameters, preserving as 
many valid models as possible, and keeping regional sizes low can 
be achieved at the Tdust /^Tdust ^ 10 level. 

A plot of Xdof distributions for 13 all-sky fixed- /3 models with 
/3 in the range 1.4 — 2.6 at 0.1 increment and Tdust /^Tdust ^ 10 
is presented in Fig.[9l The high-x^ tails of these distributions show 
that models with the largest (2.6) and smallest (1.4) values of /3 
have more fits with large Xdof • The plot of Xdof excess for these all- 



sky models is presented in Fig.|9] It shows that the all-sky /3 = 1.7 
model can fit the largest amount of data (87.6 per cent area of the 
full sky). Fig.[TO]presents sky maps of Xdof ' spatial resolution, dust 
temperature, optical depth, and the signal-to-noise of parameters 
from fitting the one-component 13 = 1.7 model. Each of the 6063 
fits presented there satisfies the Tdust /^Tdust ^ 10 requirement 
with the least amount of spatial averaging. In the parameter maps, 
if a fit has less than 10 per cent x^ probability (i.e., Xdof > 1-1-3), 
then its corresponding pixel is masked in white. Fit results for other 
13 mode ls at different levels of constraint on Tdust /^Tdust are pro- 
vided in lLiand(l201lh . 

4.3.2 /3/6/3 constraint on a free- (3 model 

We apply a similar strategy to constrain fits that use a free-/3 model. 
Instead of focusing on the signal-to-noise of the dust temperature, 
we now use signal-to-noise of the emissivity spectral index to gauge 
the amount of spectral averaging. As an example, sky maps of 
the best-fitting parameters and their signal-to-noise for the free-/3 
model with {3/513 10.0 are presented in Fig.fTT] 

Notice that with spatial averaging the dust temperature map 
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Figure 10. All-sky maps of x^of ^^^P 1^^^)' spatial resolution (top right), T^ust (second l eft) , T h i , st / (^^h u st (second right), r (bottom left), mdr/Sr (bottom 
right) of the one-component /3 = 1.7 fits that satisfy T^ust /^^^dust ^ 10.0 (from Liang 1201 ih . These maps are in Galactic coordinates MoUweide projection 
with the Galactic centre at the centre and longitude increasing to the left. The regional size map shows that majority of the fits are at the 6.7 l-D^ level. In the 
parameter maps a pixel is masked in white if it corresponds to a with less than 10 per cent probability. The group of black pixels that slant from the NEP 
to the SEP are positions where FIRAS did not provide data. 



has more consistent values at high latitudes near the Galactic poles 
than that obtained from the 7° -pixel fits in Fig. [7] The uncertainty 
of Tdust is less than 7.4 per cent as a result of the constraint on 
13/613, and the uncertainty of r has a maximum of 23.25 per cent. 
The spatial resolution of these fits are presented in an all- sky map 
in Fig. [TT] and summarized in Table [3 

The upper panel in Fig. [TJl compares Xdof distributions of 
fits that use no constraint on any parameter and those that satisfy 
P/Sp ^ 5.0, 6.7, and 10.0. The shape of the Xdof distributions 
again resemble a Gaussian, and the constrained distributions all 



centre around 0.95. This means that there is no apparent systematic 
bias in the fits and the error estimates for the data is about right. 

The plot of Xdof ^s. Galactic latitude, lower panel of Fig.[T2l 
shows that at |6| < 10°, Xdof continues to be > 1 as is the case 
of fitting all-sky one-component fixed-/3 models. At high latitudes, 
Xdof do not flare up with increasing constraint on /3/S/3, as op- 
pose to that which happens when Tdust /^Tdust requirements are 
imposed on a fixed- /3 model. 

Fig. [13] presents distributions of /3 and Tdust of one- 
component free-/3 fits that use no constraint on any parameter and 
those that satisfy 13/6/3 5.0, 6.7, and 10.0. The centres of the dis- 
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One-component free-/3 fits with (3/6(3 ^ 10.0 
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Figure 11. All-sky maps of x^of ^^^P ^^^^)' spatial resolution (top right), dust temperature (second left), emissivity spectral index (third left), optical depth 

(bottom left), and their error estimates (right column) of the best-fitting one-component free-/3 fits that satisfies (3/6^ ^ 10.0 in Galactic coordinates 

MoUweide projection with the Galactic centre at the centre and longitude increasing to the left. In the parameter maps a pixel is masked in white if it 

corresponds to a x^of ^ 1-13. The group of black pixels that slant from the NEP to the SEP are positions where FIRAS did not provide data. Notice that 

dust temperature predictions at the Galactic polar caps are noticeably lower than the rest of the sky, consistent with our knowledge that the Galactic poles have 

fewer and cooler stars. Warm regions on the temperature map correspond to known emission sources (see text). 
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Level Regional size of fits Percentage of the full sky at a regional average 







P / op ^ O.L) 


Plop ^ K). ( 


p/Op ^ iU.U 


1 


6.71 


40.01 


30.63 


22.09 


2 


60.43 


23.14 


24.53 


22.51 


3 


167.86 


11.87 


11.07 


10.60 


4 


329.00 


4.61 


8.59 


6.41 


5 


543.86 


3.53 


3.48 


5.53 


6 


812.44 


2.18 


2.78 


4.62 


7 


1134.73 


1.01 


2.12 


3.30 


8 


1510.73 




1.92 


2.59 


9 


1940.45 




1.22 


2.23 


10 


2423.88 






2.08 


11 


2961.03 






1.90 


12 


3551.89 






1.53 


13 


4196.47 






0.72 


14 


4894.76 






0.23 




Total 


87.66 


87.19 


86.34 
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Figure 8. x^of distributions (upper) and x^of Galactic latitude (lower) 
of the all-sky one-component (3 = 2.0 fits that satisfy T^ust /<^^dust ^ 
5, 10, 20 and 40, respectively. With more restrictive T^ust /<^^dust re- 
quirements, Xdof high-latitude (\b\ > 60°) fits migrate from the range 
0.7 - 1.0 to the range 0.8 - 1.3. 
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Figure 9. Upper: x^of distributions of the best-fitting one-component fixed- 
/5 models with /3 in the range 1.4 - 2.6 at 0.3 increment and satisfy 
^dust /<^^dust ^10. Lower: Percentage area of the full sky that cannot 
be fit by a one-component model with fixed emissivity spectral index (x^of 
cutoff corresponds to 10 per cent probability). This plot shows that /3 = 1.7 
models can fit the largest amount of data (87.6 per cent area of the full sky). 
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Figure 12. Upper: Distributions of x^of ^^^^ one-component free-^ fits 
that use no constraint on any parameter and those that satisfy (5/5(3 ^ 5.0, 
6.7 and 10.0, respectively. Lower: x^^f vs. Galactic latitude. 
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Figure 13. Distributions of /3 (top) and Tdust (bottom) of one-component 
free-/3 fits that use no constraint on any parameter and those that satisfy 
^/6(3 ^ 5.0, 6.7 and 10.0, respectively. Only those fits with x^of ^ l-^^ 
are included. 



tributions are at 1.80, 1.83, 1.85, and 1.88, respectively. This plot 
shows that the (3/5(3 requirement has the effect of moving (3 from 
below 1.5 to higher values. With even a moderate amount of con- 
straint on(3/3p, the range of (3 values quickly reduces to between 1 
and 3, an indication that values outside this range are rare in nature. 

Also shown in Fig. [13] the Tdust distributions do not peak at 
a single value. Instead, there is a range of most popular tempera- 
tures between 17 and 20 K. Compared to the unconstrained case, 
(3/5(3 requirement smooths out the high-temperature points and 
effectively replaces them with values at or below 20 K. 

The amounts of constraint on dust temperature, optical depth 
and emissivity spectral index for various levels of constraint on 
(3/5(3 are summarized in Table [3] In order to place adequate con- 
straint on each parameter and to keep the regional sizes low, we 
adopt (3/5(3 ^ 10.0 constraint as our standard. 



5 DISCUSSION 

5.1 Dust temperature at high Galactic latitudes 

Dust temperature is not measured directly in astronomical observa- 
tions. Instead, one derives dust temperature from intensity measure- 
ments of dust emission and a model that describes the emission pro- 
cess. That dust temperature decreases with the increase in Galactic 
latitude is supported by observations which identify thermal dust to 
be the origin of the Galaxy's diffuse emission, observations which 



show that the diffuse emission decreases from low to high Galac- 
tic latitudes following a cosecant law, observations of the cosecant 
law for dust reddening and extinction in the Galaxy with respect 
to Galactic latitude, and observations of the decreasing number of 
bright stars toward the Galactic polar caps. Substantial amount of 
work has been carried out over the last four decades, we mention 
only a few as examples: 

In reviewing g lobal properties of dust in the Galaxy, 
IStein & Soifej ([1983; point out that "the far-infrared luminosity in 
the Galaxy is associated with regions where stars are now or were 
recently forming," and so "even in the galactic center region, much 
of the far-infrared luminosity is associated with ongoing or recent 
star formation activity." Using far-infrared luminosity as a measure 
of the current rate of star formation, dust acts as "a frequency con- 
verter which absorbs the short wavelength photons emitted by the 
newly formed stars and re-emitting the energy at far-infrared wave- 
lengths." 

That dust grains are heated by the general interstellar radi- 
ation field (ISRF) which follows a cosecant law to the Galac- 
tic poles and has a stellar origin have been confirmed by many 
research groups. For instan ce, by fitting s pectra of diffuse emis- 
sion from the Galactic plane iMezger et all (Il982h and lMathis et al.l 
(1983) identify two main contributors to the diffuse emission at 
wavelengths ^ 20 /xm: dust heated by O stars (60%) and dust 
associated with diffuse atomic intercloud gas heated by the gen- 
eral interstellar radiation field (40%). They conclude that dust as- 
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Table 3. Parameters and their uncertainties of the one-component free-/3 fits with x^^f ^1.13 
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^dust 


<^^dust 






(K) 


(K) 




min 


max 




^ 5.0 


10.11 


23.83 


2.67 


^ 6.7 


10.12 


22.69 


1.95 


^ 10.0 


13.69 


22.69 


1.26 




Constraint on (5/5(5 










min 


max 




^ 5.0 


1.08 


4.71 


0.91 


^ 6.7 


1.08 


4.80 


0.71 


^ 10.0 


1.24 


3.13 


0.31 




Constraint on (5/5(5 




r 


5t/t 




xlO-^ 


(per cent) 




min 


max 




^ 5.0 


0.33 


46.15 


62.83 


^ 6.7 


0.45 


46.15 


49.02 


^ 10.0 


0.61 


46.15 


23.25 



sociated with clouds containing no luminous sources of heating 
such as OB stars contributes less than 7% of the total diffuse far- 
infrared/sub-millimeter emission. In 1986, Cox, Kriiegel & Mezger 
showed that the total luminosity of the re-radiated dust emission 
from the Galactic disk is about 40% of the total stellar luminos- 
ity. Their best-fitting dust model predicts that ^ 90% of the total 
Galactic infrared luminosity come from three components: a cold 
dust (T ^ 15 — 25 K) component associated with atomic hy- 
drogen, a very cold dust (< T >^ 14 K) component associated 
with molecular hydrogen and located inside quiescent molecular 
clouds, and a warm dust (T ^ 30 — 40 K) component associated 
with ionized gas in extended low density Hii regions heated by O 
and B stars. The first two components are heated by the general 
ISRF generated by both young and old stellar populations. They 
dominate the sub-millimeter part of the spectrum and have a to- 
tal luminosity of ^ 37% of the total Galactic infrared luminosity. 
The luminosity of the warm dust component, on the other hand, ac- 
counts for ^ 50% of the total Galactic IR luminosity. Using the 
first all-sky infrared measurements taken by the Infrared Astro- 
nomical Satellite (IR AS, IRAS Explanatory Supplement (1988)), 
iBoulanger & Peraultl ( 119881) show that Galactic emission exists at 
Galactic latitude \b\ > 10° at 12, 25, 60 and 100 /xm and that the 
emission comes mainly from dust heated by the ISRF associated 
with atomic gas and diffuse ionized gas. More specifically, they 
show that the Galactic latitude profiles of the IRAS all- sky mea- 
surements at 60 and 100 /xm follow a cosecant law from 10° to the 
poles. This means that dust emission decreases from low to high 



latitudes in the manner as would be expected for a plane-parallel 
configuration. 

Not only does the strength of the ISRF decreases with in- 
creasing latitude, but the amount of dust at high Galactic latitudes 
also decreases following a cosecant law. Distribution of dust at 
high Galactic latitudes can be infer red from studies of dust red - 
dening and extinction. For example, iMcClure & CrawfordI (Il97lh 
use 207 late-type stars (K giants) near the north Galactic pole to 
calculate a reddening law and find that it is consistent with the 
cosecant law (0.06 • esc 6), where the Sun is assumed to be in a 
cylindrical hole of an absorbing layer in the Galactic plane. They 
conclu de that there is almost no abs orption from the pole down 
to 50°. lde Vaucouleurs & Butal (Il983h studies high-latitude extinc- 
tion in extragalactic objects using five independent methods: counts 
of faint galaxies, counts of distant galaxy clusters, mean surface 
brightness of bright galaxies, hydrogen-luminosity ratio measured 
by the Hi index for bright galaxies, and total color indices of bright 
galaxies. They show that extinction derived from each of these 
methods supports the cosecant law from \b\ > 20° to the poles 
in each Galactic hemisphere. 

Finally, we know that stars are the power houses in the Galaxy 
for both the ISRF and interstellar dust grains and that all- sky sur- 
veys have shown that the Galactic p olar regi ons have fewer and 
cooler stars (sky maps in IRAS Expl a nator^SuDDlemeni Il988l : 
lArendt et all Il998l : IVoges et all 1 19991 : iGudennavar et all 120121) . 
With the observations of less intense emission, fewer dust parti- 
cles and fewer sources of energy at the Galactic polar regions, we 



© 2011 RAS, MNRAS 000,[TH20l 



14 Z Liang, D. J. Fixsen and B. Gold 
301: 




TOO 



GaEoctic lotEtude (*) 



Figure 1 5. Thi 1st vs. Galact ic latitude for the two -component (3 = 2.0 
model of iReach et alJ Il995h (black horizontal bars), the two-c omponent 
(3i = 2.70 and (32 = 1.67 Model #8 ofjRnkbeiner etaD (ll999h (purple), 
the one-component 13 = 1.7 model of iLiana C2011«) (red), and the one- 
component free- (3 model from this work (blue). 



30 
25 
20 
15 
10 
5 




2.0, lbl>60° 
, lbl<50° 
1.7, lbl>60° 
, lbl<60° 

□ this work, |8/(5j8 > 10.0, lbl>60° 
, lbl<60° 



□ Reach et al. (1995), 

X 

C Liang (201 1), / 



10 



Figure 16. T, 



100 
T/(5T 



1000 



of iReach et al 



; t VS. T h 1 1 St /<^?ri 11 St, for the two-component ^ = 2.0 model 
(Il995h (green and black), the one-component (3 = 1.7 
model (orange and light blue), and the one-component free-/3 model (ma- 
genta and blue). Squares represent fits with |6| > 60° , and crosses represent 
fits with |6| < 60°. 



thus expect that dust temperature in those regions to be lower ac- 
cordingly. 

5.2 Comparing dust temperature predictions of different 
models 

The crucial test to any model is comparison with observations. In 
this section, we present all- sky maps of dust temperature obtained 
from the literature and those from our study. We argue that since 
dust temperature distribution should generally decrease with the 
increase in Galactic latitude, models that fail such a test are un- 
physical. 

Fig. [14] presents dust tempe rature predic ti ons b y the two- 
component 13 — 2.0 model of iReach et aP (Il995h . the two- 
component 01 = 2.70 and (32 = 1.67 Model #8 of 
iFinkbeiner etaP (Il999h . the one-component f3 = 1.7 model with 
Tdust/^^dust ^ 10.0, and the one-component free-/3 model with 
^ 10.0. The last two models are the ones derived in Section 
|4] We do not h ave Planck Collaboration's dust temperature predic- 
tions JPlanck Collabora tion 2011a), but since their one-component 
dust model has the emissivity spectral index at 1.8, we use the tem- 
perature map of our one-component /3 = 1.7 model to approximate 
their results. 

Notice the difference in dust temperature predictions at the 
Galactic polar caps relative to lower latitude regions in these maps. 
In the first three cases, where emissivity spectral indices of the 
models are fixed, temperature estimates around the Galactic polar 
caps are noticeably higher or at about the same level compared to 
those at lower latitudes. On the other hand, the free-/3 model pre- 
dicts the polar caps to have the lowest temperature of the entire 
Galaxy. 

To perform the comparison in a more quantitative manner, we 
plot Tdust vs. Galactic latitude for the different models in Fig. 
[T5l One may notice that data points corresponding to the one- 
component /3 = 1.7 model flare up at |6| > 30° most dramatically 
among all models. One may also notice that the two temperature 
components of Model #8 share the same distributions and that their 
values are approximately independent of Galactic latitude. In the 
southern hemisphere, one may notice that the two dust components 
of the Reach model consistently have higher temperatures at the 
Galactic polar caps than at mid latitudes. In fact, temperature of 



the cold dust at the southern polar cap is the highest that cold dust 
experiences throughout the entire Galaxy. In comparison to these 
three models, the free-/3 model is the only model which predicts a 
decrease in dust temperature with respect to Galactic latitude. 

We investigate the possibility that dust temperature predic- 
tions are biased positively when signal-to-noise of data decreases. 
Fig. US] plots Td ust as a function of Tdust /<^^dus t as gi ven by the 
two-component = 2.0 model of iReach et al ] (Il995h . the one- 
component 15 = 1.7 model, and the one-component free-/3 model. 
Model #8 is not included because FDS study does not provide er- 
ror estimates for their models. Notice that while signal-to-noise 
of the fits vary by over one, sometimes two, order of magnitude, 
Tdust values change by less than 70 per cent, so Tdust is approxi- 
mately independent of the signal-to-noise of the fits. More specifi- 
cally, within the set of points for each model, in comparison to the 
(Tdust|^l>gQO, Tdust/(^Tdust|^|>gQo) pairs, the existence of lower 

Tdust /^Tdust and lower Tdust pairs or higher Tdust /^Tdust and 
higher Tdust pairs shows that the data do not support a definitive 
causal relation between low signal-to-noise and high temperature 
prediction. Based on these observations, we conclude that any po- 
tential bias to the fits due to low signal-to-noise cannot be the main 
reason for high temperature predictions seen at the Galactic polar 
caps by fixed- /3 models. 

Since we have checked temperature predictions of models 
with sp ectral index a t the range of most likely values between 1 .4 
and 2.6 (lLiandl20T]]) . and they all exhibit higher temperature at the 
Galactic polar caps than at lower latitudes, we suspect that such a 
phenomenon is an inherent artifact resulted from fixing the emis- 
sivity spectral index at a particular value. As a result, we argue that 
the temperature predictions of these fixed- /3 models are unphysical. 

In addition to predicting reasonable temperature distribu- 
tion at the Galactic polar caps, the free-/3 model predicts re- 
gions with known heating sources to have higher temperature than 
their surrounding environment. For instance, we can identify, off 
the Galactic plane, high temperature regions that correspond to 
the Ophiuchus region, {3 Centauri, the Large Magellanic Cloud 
and the Orion region. At the Galactic plane, model predictions 
show that dust temperature decreases with increasing Galactocen- 
tric di stance, consistent wi th observations of the DIRBE experi- 
ment (ISodroski et al.lll99"4l) . On top of this temperature gradient 
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Figure 14. Dust temperature maps by the two-component /3 = 2.0 model of iReach et alJ (Il995h (first row), the two-compo nent (3 i =2.70 and (32 = 1.67 
Model #8 of Finkbeiner et al. (1999) (second row), the one-component (3 = 1.7 model with T^ust/^^dust ^ 10 of Lian^ 1201 ih (third row), and the one- 
component free-/3 model with (3/513 ^ 10 of this work (last row). Data of the Reach model were obtained from their manuscript; the Galactic plane of both 
maps have been masked due to incomplete data for the region. Data of the Finkbeiner model were generated from their IDL code "predict_thermal.pro." Pixels 
in both the /3 = 1.7 and the free-/5 maps are masked in white if they correspond to a x\of > 1.13. 
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Figure 17. Averaged dust spec tra (orange) and pr edictions of the two- 
component ^ = 2.0 m odel of lReach et aP (Il995h (blue), Model #8 of 
i Finkbeiner et al.l (Il999h (magenta), one-component /3 = 1.7 model of 
iLiand ( l201lh (green), and one-component free-/3 model (red). 



Figure 18. Difference between the averaged 214-channel dust spectra and 
model predictions. Color-model assignment is the same as in Fig.fTT] The 
differences are obtained by first subtracting model predictions from each 
7° dust spectra, then dividing the results by measurement uncertainties, and 
finally averaging the differences. 



the temperature map shows the Outer Galaxy to be asymmetrically 
interspersed with strong infrared emission regions, such as NGC 
7538 at / = 111°, Cas A at 111?7 < / < 112? 1, W3/4/5 at 
132? 5 < / < 138? 5, Cygnus region at 75° < / < 90°, and the 
Gum Nebula at 187? 5 < / < 193? 5. These regions at the Outer 
Galaxy are currently masked under the criterion that Xdof of 
fit be < 1.13. For these regions we believe that complex emission 
sources in addition to thermal dust have contributed to our prepared 
spectra. This challenges our model's assumption that only thermal 
dust can contribute to the continuous emission. If true, then our 
model is no longer sufficient. 

Fig. [it] present samples of fits to the average dust spectra by 
each of the aforementioned models at different Galactic latitudes. 



While most models are able to simulate the data with various de- 
grees of success, the free-/3 model is consistently able to trace the 
data points closely at all Galactic latitudes. Fig. [Tsl compares the 
normalized difference between data and intensity predictions for 
the same set of models and the same sky regions. It shows that the 
free- 13 model is able to match the data within 1-sigma for majority 
of the channels while the other models are less consistent at pre- 
dicting close to measurement values across frequency and Galactic 
latitude. That the normalized differences of the /3 = 1.7 model are 
comparable to those of the free-/3 model demonstrates the limita- 
tion of relying on a single number to evaluate a multifaceted opera- 
tion, such as using the normalized difference or the to differen- 
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tiate models. Here, the free-/3 model and the {3 — 1.7 model have 
about the same value, but one of them has a strong connection 
with other facts we know of the Galaxy and the other does not. 

5.3 Relation between dust temperature and emissivity 
spectral index 

The existence of an anti-correlation between the emissivity spectral 
index and dust temp erature has been a subje ct of debate in the last 
decade. For example. [Veneziani"et did not find a clear in- 

verse correlation between the two p arameters in BOOMERANG 's 
measurements of the Galactic cirrus. IPlanck Collaboration! (l2011ah 
derived all- sky dust temperature using a fixed value of the emis- 
sivity spectral index due to degen eracy between ten iperature 
and spectral index i n a fre e-/3 fit dShettv et all l2009al lbl). And 
IPlanck Collaboration! (l2011ch carried out analysis of dust proper- 
ties at high Galactic latitudes using the same fixed spectral index 
model in order to maintai n consistency within the C o llabor ation. 

On the other hand, IPupac et all (1200 iL l200l l2003h of the 
ProNaOS experiment reported evidence of the inv erse correla- 
tion in their observations of M42, Ml 7 and NGC89L lbesert et all 
(l2008h of the Archeops experiment confirmed the correlation on 
300 Galactic point sources, ■FaradisUioToh confirmed the existence 
of the correlation in two Hi-G AL fields at the Galactic plane, and 
iFlanck Collaboration! feOllbl) concluded that observation of ob- 
jects in the Early Cold Core Catalogue were not consistent with 
a constant value of the spectral index over all temperature. 

The strongest argument for the existence of a temperature 
dependent power-law emissivity spectral in dex comes from lab- 
oratory measurements. lAdadze et all ( Il996h measured absorption 
spectra of crystalline and amorphous silicate grains at 0.7 - 2.9 mm 
at 1.2 - 30 K and found that the wavelength dependent absorption 
coefficients of the amorphous grains can be modeled with power 
laws, and the exponents of which vary with temperature. Agladze et 
al. attributed the correlation to a two-level population eff ect found 
in the low-lying tunneling states in bulk glasses. Menn ella et al.l 
(Il998h measured the absorption coefficients of crystalline and 
amorphous carbon and silicate grains at 20 /xm - 2 mm at 4 - 295 K 
and found that the absorption coefficient followed a power law with 
the spectral index depending on grain temperature. They attributed 
the observed pattern to the dominance of two-phonon difference 
processes. In 2001, Chihara et al. measured far-infrared absorption 
spectra of crystalline circumstellar dust analogues between 20 - 
80 fim at 4.2 K and at 295 K. They also reported a higher value 
of the spectral index at 4.2 K than at 295 K. In 2005, Boudet et 
al. measured the mass absorption coefficients of amorphous sili- 
cates between 100 fim - 2 mm at 10 - 30 K and found that the 
power law spectral index had a strong temperature and wavelength 
dependence. More specifically, the anti-correlation between spec- 
tral index and grain temperature was more pronounced at 500 fim 
- 1 mm than at 100 - 250 /xm. Comparing their measurements 
with the temperature- spectral index relation deduced from synthe- 
sis of the FroNaSO observations, Boudet et al. found a comparable 
anti-correlation between temperature and spectral index between 
observations and laboratory measurements when the spectral in- 
dex is between 1 and 2.7. They attributed the origin of the relation 
to the tunneling processes in two level systems. Bringing together 
evide nce from laboratory measurements and astronomical observa- 
tions jMenve t al. (2007) proposed a model that relates the emission 
spectrum's dependence on temperature and internal structure of the 
grains. Additional evidence in support of the explanation that in- 
ternal grain structures give rise to the apparent correlation between 
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dust te mperature and spectral index came in when lCoupeaud et aP 

l|) reported that mass absorption coefficients of synthesized 
analogues of interstellar amorphous silicate grains decreased with 
temperature and that the local spectral index anti-correlated with 
grain temperature. These laboratory results strengthen our confi- 
dence that the anti-correlation between dust temperature and spec- 
tral index seen in astronomical observations is not spurious but in- 
stead grounded on the physical properties of interstellar dust grains. 

To test the possibility of an inverse correlation between emis- 
sivity spectral index and dust temperature in our data, we fit the 
best-fitting /3 and Tdust values of the free-/3 model with 13/ 6/3 ^ 
10.0 to the hyperbolic function (3 = 1/(6 -\- u ■ Tdust) and obtain 
S = -0.510 ± 0.011 and uj = 0.059 db 0.001 with a Xdof = 0-99. 
Fig. [19] presents the data and the Tdust-/^ relations identified by 
this and other works. Notice that our best-fitting Tdust -/^ relation 
overlaps with all but that of the Hi-GAL field at / = 30°. Since 
the (Tdust, 13) pairs we used to constrain the hyperbolic function 
come from all over the sky except for about 13 per cent of the total 
sky area in the Inner Galaxy, the overlap with other experiments' 
results confirms that our Tdust -/^ correlation has a broad represen- 
ta tion. That it is w ithin one-sigma of the Tdust- f3 relation identified 
m |Paradisl ( [20TQh for the / = 30° field may have to do with one or 
both of the following reasons. First, the physical properties of dust 
may have changed from the Inner Galaxy to regions far away from 
active star formation. If so, then our set of (Tdust, 13) values is not a 
good representation of dust properties in the / = 30° field since fits 
of the Inner Galactic spectra are not included. Secondly, one may 
understand the isolation o f the curve for Hi-GAL field at / = 30° 
from explanation given in IParadisI tOld) , which points to the fact 
that data from the / = 30° field challenge their model's single- 
temperature assumption. Finally, that our set of (Tdust, 13) values 
marginally overlap with the relation from the Plan ck study of cold 
cores in the Galaxy dPlanck Collaboration! |201 Ibl) may have to do 
with the fact that there is only a smal l overlap between the two data 
sets: While the 7608 cold cores in the lPlanck Collaboration! (!2011bl) 
analysis largely distribute along the Galactic plane, our model's 
target regions are at mid and high latitudes and part of the Outer 
Galaxy. In addition, the 7608 cold cores form a selective group of 
objects. To be selected, an object has to be a compact source colder 
than the surrounding envelope and a 17-K diffuse Galactic back- 
ground. This last criterion would have excluded over 73 per cent 
of our data points. Fig. |20] presents the constant contours, cor- 
responding to confidence intervals with 68, 95 and 99.99 per cent 
probability. It shows that a temperature independent (3 model is re- 
jected by the data with more than 99.99 per cent probability, an 
indication for a real dependence between Tdust and (3. 

We present a summary of the best-fitting 6 and cj for other 
levels of constraint on the fits in Table |4] The exact functional form 
of the temperature- dependent spectral index cannot be established 
by our data alone. However, that a physically sensible fit to the all- 
sky data requires a free-/3 model and that an inverse relation fits the 
13 and Tdust values well give strong support to the anti-correlation 
between the emissivity spectral index and dust temperature. 

The observed correlation between emissivity spectral index 
and dust temperature in our data is likely to be intrinsic rather 
than from noise. In 2009, Shetty et al. showed that the accuracy 
of parameter estimation depends on the wavelength range of the 
spectra used in model fitting. They show that for dust temperature 
^ 20 K, fits using spectra at 100 — 600 /xm are not sensitive to 
noise. However, if 40 K ^ Tdust ^ 80 K, a 5 per cent noise in 
a data set of different Tdust and the same [3 may generate a pat- 
tern in the best-fitting {3 vs. Tdust plot that can be confused with a 
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Table 4. Best-fitting parameters to the relation (5 = l/{5 -\- uo • T^ust) between emissivity spectral index and dust temperature 



constraint on fits 


5 




number of data points 


^dof 


no constraint 


-0.305 ±0.014 


0.048 ± 0.001 


5421 


0.72 


(3/5(3 ^ 5.0 


-0.399 ±0.013 


0.053 ± 0.001 


5386 


0.83 


^/5(3 ^ 6.7 


-0.441 ±0.012 


0.055 ± 0.001 


5357 


0.88 


(3/5^ ^ 10.0 


-0.510 ±0.011 


0.059 ±0.001 


5305 


0.99 




Figure 19. Data and the best-fitting /3 — T^ust relations from this and other 
observations. The black crosses and error bars represent best- fitting values 
from one-component free-/5 fits that satisfy (3/5(3 ^ 10.0. Both data and 
errors are used in fitting the hyperbolic function (3 = 1/(5 -\- u • T^ust)- 
Only fits with x^of ^ ^-^^ included. Also presented are the best- 



fitting /3 — Tdust relation of the ProNaOS experiment (blue, DupacetajJ 
l2003h . the Archeops e xperimen t (magenta, Desert et al. 2008), the two 
Hi-GAL fields (green, Paradis 2010), and the Planck mission (orange, 
IPlanck Coll aboration 2011b). Dashed curves represent 1-sigma uncertain- 
ties of the respective models. 
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Figure 20. Constant contours corresponding to confidence intervals with 
68, 95 and 99.99 per cent probability. We fit the relation /3 = 1/(5 ± • 
^dust ) to the best-fitting T^ust and /3 values from the one-component free- 
ze model. The white cross represents location of the minimum • 



true anti-correlation in the two parameters. Since our spectra, com- 
ing at 100 — 4400 /im, well sample the Rayleigh- Jeans tail of the 
peak of the ^ 20 K modified blackbody spectrum, our best-fitting 
parameters are good approximates of the true values. 

As the quality of data improves, the free- (3 model will allow 
future studies to pinpoint the values of (3 and Tdust with greater ac- 



Figure 21. Scatter plot of /3 /5^ and T^ust /^^dust • The blue, green, yellow 
and red data points represent values of and T^ust /(^^dust from fits 
that use no constraint on any parameter and those that satisfy (3/5(3 ^ 5.0, 
6.7 and 10.0, respectively. Only fits with X^^f ^1.13 are included. 



curacy. In Fig.|2T]we present (3/6(3 vs. Tdust /^^dust of the all-sky 
collections of fits with different amounts of constraint on (3/6(3. li 
shows that the uncertainty in (3 and Tdust are directly proportional 
while (3 is more susceptible to errors than Tdust - That is, with more 
sensitive data one will be able to reduce uncertainty in (3 and Tdust 
at the same time. 



6 CONCLUSIONS 

We have examined models of interstellar dust thermal emission at 
far-infrared and millimeter wavelengths. Starting with deducing a 
new and improved version of dust spectra from FIRAS calibrated 
measurements and unifying them with measurements of the DIRBE 
and the WMAP, one-component dust models with fixed and vari- 
able emissivity spectral index are fit to 6063 214-channel spectra at 
fixed and variable spatial resolutions. 

We show that the free- (3 model can predict more physically 
motivated dust temperature than models with fixed emissivity spec- 
tral index. With observations showing that the intensity of dust 
emission and the amount of dust decrease with the increase in 
Galactic latitude and that fewer and weaker stars exist at higher 
Galactic latitudes, dust temperature is expected to be lower at the 
Galactic polar caps than at lower latitudes. Among all- sky dust 
models found in the literature and the one-component dust models 
constructed in our study, we demonstrate that only the free- (3 model 
passes this test. All fixed- /3 models predict higher or the same tem- 
peratures at the Galactic polar caps as at lower latitudes. 

Our best dust model is the one-component free- (3 model with 
(3/6(3 10.0. It fits dust spectra over 86 per cent area of the full 
sky. Dust temperature is predicted to be 13.69 - 22.69 (±1.26) K, 
emissivity spectral index to be 1.2 — 3.1 (±0.31), and optical depth 
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to range 0.61 - 46.15 x 10~^ at = 900 GHz (Ao = 333/im) 
with a 23.25 per cent uncertainty. Our model has the 7° angular 
resolution of FIRAS and uses variable spatial averaging at high 
latitudes. It supports the interpretation of an anti-correlation be- 
tween emissivity spectral index and dust temperature. We fit the 
relation P = + uj ■ Tdust) to the /3 and T values and obtained 
5 = -0.510 ± 0.011 and cj = 0.059 =b 0.001. A temperature in- 
dependent 13 model (cj = 0) is rejected by the data with more than 
99.99 per cent probability, an indication that a real dependence be- 
tween Tdust and {3. 

The free-/3 model can be used to remove dust contamination 
in CMB measurements. It can also serve as an all- sky extended- 
frequency-range reference for other experiments. With higher res- 
olution data sets, such as those from Planck, one will be able to 
further refine our knowledge of interstellar dust grains and their 
interaction with other constituents in the Galaxy. 
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